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Abstract.  The  problem  of  heat  transfer  and  temperature  distribution  in  a  binary  mixture  of  rarefied  gases 
between  two  parallel  plates  with  different  temperatures  is  investigated  on  the  basis  of  kinetic  theory.  Under 
the  assumption  that  the  gas  molecules  are  hard  spheres  and  undergo  diffuse  reflection  on  the  plates,  the 
Boltzmann  equation  is  analyzed  numerically  by  means  of  an  accurate  finite-difference  method,  in  which 
the  complicated  nonlinear  collision  integrals  are  computed  efficiently  by  the  deterministic  numerical  kernel 
method.  As  a  result,  the  overall  quantities  (the  heat  flow  in  the  mixture,  etc.)  as  well  as  the  profiles  of 
the  macroscopic  quantities  (the  molecular  number  densities  of  the  individual  components,  the  temperature 
of  the  total  mixture,  etc.)  are  obtained  accurately  for  a  wide  range  of  the  Knudsen  number.  At  the  same 
time,  the  behavior  of  the  velocity  distribution  function  is  clarified  with  high  accuracy. 


INTRODUCTION 

The  problem  of  heat  transfer  and  temperature  distribution  in  a  rarefied  gas  between  two  parallel  plates  with 
different  temperatures  is  one  of  the  classical  problems  in  rarefied  gas  dynamics,  and  a  large  number  of  theoretical 
and  experimental  works  have  been  devoted  to  this  problem,  especially  in  the  case  of  single-component  gases 
(see,  e.g.,  [1-4]  and  the  references  cited  in  [3,4]).  Early  theoretical  works  covering  a  wide  range  of  the  Knudsen 
number  were  mainly  based  upon  either  moment  and  variational  methods,  containing  arbitrary  assumptions  on 
the  form  of  the  velocity  distribution  function,  or  numerical  analysis  using  model  Boltzmann  equations.  Only 
in  1989,  Ohwada  et  al.  [3]  reported  an  accurate  numerical  solution  of  the  linearized  Boltzmann  equation  for 
a  hard-sphere  gas  in  the  case  of  a  small  temperature  difference  between  the  plates.  Their  solution  method 
was  a  finite-difference  method,  in  which  the  collision  integral  was  computed  efficiently  as  well  as  accurately  by 
the  numerical  kernel  method  developed  by  Sone  et  al.  [5] .  Subsequently,  Ohwada  extended  the  method  to  the 
collision  integral  of  the  full  Boltzmann  equation  in  his  shock-structure  analysis  [6]  and  then  applied  it  to  the 
heat-transfer  problem  for  a  nonsmall  temperature  difference  between  the  plates  [4,7]. 

As  for  the  case  of  binary  gas  mixtures,  the  accumulation  of  the  results  is  not  satisfactory,  though  some 
analyses  (by  means  of  a  moment  method)  as  well  as  experiments  were  performed  in  1970’s  [8,9].  In  the  present 
study,  therefore,  we  investigate  the  heat-transfer  problem  for  a  binary  mixture  of  hard-sphere  gases  on  the  basis 
of  the  full  Boltzmann  equation  for  a  large  temperature  difference,  aiming  to  provide  an  accurate  numerical 
solution  that  can  be  regarded  as  a  standard  for  the  problem.  Recently,  the  present  authors  have  extended 
Ohwada’s  numerical  kernel  method  for  the  nonlinear  collision  integral  to  the  case  of  binary  mixtures  in  their 
study  of  shock  wave  structure  [10].  The  same  method  is  employed  in  the  present  analysis. 


PROBLEM 

Consider  a  rarefied  mixture  of  two  gases,  say  components  A  and  B,  in  the  domain  0  <  Xi  <  D  between  two 
parallel  plane  walls  at  rest,  where  Xi  is  a  rectangular  coordinate  system  in  space.  Let  the  wall  at  Xi  =  0  be 
kept  at  temperature  T/  and  that  at  Xi  =  D  at  temperature  T//.  Investigate  the  steady  behavior  of  the  mixture 
(temperature  distribution,  heat  flow,  etc.)  on  the  basis  of  kinetic  theory  under  the  following  assumptions: 
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(i)  The  molecules  of  each  component  axe  hard  spheres,  and  the  interaction  between  two  gaseous  molecules  is 
the  complete  elastic  collision. 

(ii)  The  molecules  of  each  component  are  reflected  according  to  the  diffuse  reflection  condition  on  the  walls. 


BASIC  EQUATION 


Let^  =  (Cii  C2,  Ca)  be  the  molecular  velocity  and  F"(Xi,  4)  the  velocity  distribution  function  of  the  molecules 
of  a-component  {a  =  A,  B).  The  Boltzmann  equation  in  the  present  problem  is  written  as 


E  («  =  AS),  (1) 

^  0=A,B 

=  [(4“)V2]  I mengien  -  HUgim^  •  (2) 

4“  =  +  di)/2,  +  m^).  (4) 


Here,  m"  and  are  the  mass  and  diameter  of  a  molecule  of  a-component;  is  the  integration  variable  for 
e  is  a  unit  vector,  =  c!^*id$+2d^*3,  and  dQ  is  the  solid-angle  element  around  e;  the  domain  of  integration 
is  the  whole  space  of  4*  and  all  directions  of  e. 

The  boundary  condition  on  the  walls  (Xi  =  0  and  D)  is  expressed  as  follows:  For  £  •  n  >  0, 


V  Klyj  /  J^^.nKO 


(5) 


where 


Tw=Ti,  n  =  (1,0,0),  at  Xi  =  0,  and 


Tyj  =  Tii,  n=(— 1,0,0),  at  Xi  =  U, 


(6) 


and  k  is  the  Boltzmann  constant. 

If  we  rewrite  the  equation  and  boundary  condition  in  a  dimensionless  form,  we  find  that  the  problem  is 
characterized  by  the  following  five  parameters:  jm^,  d^/d^,  Tu/Tj,  and  Kn.  Here,  n“^  is  the 

average  molecular  number  density  of  a-component  in  the  domain  0  <  Xi  <  D,  and  Kn=  Iq/D  is  the  Knudsen 
number,  where  Iq  =  [\/27r (d:^)^nav]^^  is  the  mean  free  path  of  the  molecules  of  H-component  when  it  is  in  the 
equilibrium  state  at  rest  with  number  density  Uav  = 


NUMERICAL  ANALYSIS 


We  first  note  that  in  the  present  problem  we  can  seek  the  solution  in  the  form  F“(Xi, ^1, 77),  where  rj  = 


(€1 


^1)1/2  We  analyze  Eqs.  (l)-(6)  numerically  by  means  of  an  iterative  finite-difference  method.  The  key 


TABLE  1.  Heat  flow  qt  =  {qi,  0,  0)  of  the  total  mixture  for  Tn/Tj  =  2.  Here, 


po  =  kUavTi  is  the  reference  pressure. 


^av  l^av 

Kn 

vnP  jm^  =  0.5, 

=  1 

jm^  =  0.25,  d^/d^ 

=  0.5 

gi/[po(2fcTj/m'^)^^^] 

A  (%) 

qi/\po{2kTi 

A  (%) 

0.1 

0.1 

-0.184 

0.45 

-0.207 

0.72 

0.1 

1 

-0.509 

0.19 

-0.547 

0.19 

0.1 

10 

-0.656 

0.049 

-0.693 

0.047 

1 

0.1 

-0.209 

0.34 

-0.370 

0.67 

1 

1 

-0.589 

0.15 

-0.814 

0.13 

1 

10 

-0.763 

0.047 

-0.966 

0.036 

10 

0.1 

-0.245 

0.34 

-0.659 

0.19 

10 

1 

-0.677 

0.10 

-1.124 

0.075 

10 

10 

-0.871 

0.038 

-1.244 

0.014 
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issue  in  the  analysis  is  an  accurate  and  efficient  conaputation  of  the  complicated  collision  integral  using 
the  discrete  values  of  at  the  grid  points  )  in  the  (Xi , , ??)  space.  For  this  purpose, 

we  expand  F“  at  Xi  =  as 

&  =  ii 

hi 

where  is  a  localized  basis  function  that  is  sectionally  quadratic,  takes  unity  at 

and  is  nonzero  only  in  its  neighborhood;  Li{y)  is  the  Laguerre  polynomial  in  y  of  order  1.  The  coefficients  a“/*^ 
are  determined  in  such  a  way  that  Eq.  (7)  coincides  with  at  the  grid  point  If  we  substitute 

Eq.  (7)  into  the  collision  integral  it  is  expressed  as  a  linear  combination  of  the  collision  integrals 

for  the  functions  of  the  form  (77“)^™$“.  The  latter  collision  integrals  are  independent  of  and  therefore  can 
be  computed  beforehand  (numerical  collision  kernel).  Once  the  numerical  kernel  is  prepared,  the  computation 
of  the  collision  integral  in  each  iteration  step  is  reduced  to  simple  products  and  sums  of  matrices.  In  this  way, 
high  efficiency  in  the  computation  of  the  collision  integral  is  attained  (see  [6,10]  for  the  details). 


F“(xW,a,^)  =  exp(-^) 


RESULT  OF  ANALYSIS 

The  computation  was  carried  out  for  T///T/  =  2,  (m®/rn^,  d^/d^)  =  (0.25, 0.5)  and  (0.5, 1),  =  0.1, 

1,  and  10,  and  Kn=  0.1,  1,  and  10.  To  show  the  result,  we  denote  by  n“  the  molecular  number  density 
of  a-component  (a  =  A,B)  and  by  T  and  q,  =  (q'i,0,0)  the  temperature  and  the  heat  flow  of  the  total 
mixture,  respectively  [i.e.,  n“=/ {3knT,  2qi)=f(l,  +  m^F^)d^,  where  n  = 

and  d^=d^id^2dCa].  Note  that  the  flow  velocity  of  each  component  vanishes  identically  and  the  heat  flow  qi  is 
independent  of  Xi  in  the  present  problem  because  of  the  conservation  of  mass  and  that  of  energy. 

The  values  of  qi  in  all  the  cases  are  shown  in  Table  1,  where  po  =  kuavTi  is  a  reference  pressure.  The 
numerical  result  of  qi  varies  slightly  with  Xi  because  of  numerical  error.  Its  average,  say  qiavt  over  0  < 
Xi  <  D  is  shown  as  qi  in  the  table.  The  maximum  variation  of  qi  over  0  <  Xi  <  D  relative  to  qiav- 
A  =  maxjq'i  —  i7ia^;|/|qiai;|5  which  gives  a  good  measure  of  accuracy  of  the  computation,  is  shown  in  percentage 
in  Table  1.  Figure  1-3  show  the  profiles  of  the  number  densities  and  and  of  the  temperature  T  for 
the  case  of  /m^  =  0.25,  d^/d^  =  0.5:  Fig.  1  is  for  =  0.1,  Fig.  2  for  n^y/n^y  =  1,  and  Fig.  3  for 

n^vln^y  =  10.  The  smaller  molecules  (the  molecules  of  S-component)  have  a  larger  mean  free  path.  Since  Kn 
is  based  on  the  average  number  density  of  the  total  mixture  and  on  the  diameter  of  the  larger  molecules  (the 
molecules  of  A-component),  the  effective  Knudsen  number  at  the  same  Kn  is  larger  for  larger  values  of  n^y/n^y 
Therefore,  the  temperature  jump  on  the  walls  at  the  same  Kn  is  larger  for  larger  n^y/n^y  In  Figs.  1-3,  the 


FIGURE  1.  Profiles  of  the  number  densities  and 
mF IvnA  =  0.25,  dF/d^  =  0.5,  and  riav/nav  =  0.1.  Here, 
method,  and  •  and  o  that  by  the  DSMC  method. 


and  of  the  temperature  of  the  total  mixture  T  for 
—  and  —  indicate  the  result  by  the  finite-difference 
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corresponding  result  obtained  by  the  direct  simulation  Monte  Carlo  (DSMC)  method  [11]  is  also  shown  for 
comparison.  The  DSMC  result  shows  good  agreement  with  the  finite-difference  result. 

The  velocity  distribution  functions  at  two  points  near  the  walls  in  the  cases  corresponding  to  Figs.  1-3 
are  shown  in  Figs.  4-6,  respectively.  That  is,  F°^  at  'q{2kTi =  0.15  and  1.35  are  shown  as  functions  of 
^1.  In  the  case  of  free- molecular  gas  (Kn  =  oo),  the  velocity  distribution  functions  for  any  Xi  are  discontinuous 
at  =  0.  For  large  Kn  (Kn  =  10),  though  the  discontinuity  vanishes  because  of  the  molecular  collision,  the 
gradient  near  =  0  is  still  very  steep.  The  change  around  =  0  becomes  milder  as  the  Knudsen  number 
decreases.  The  corresponding  result  by  the  DSMC  method  is  also  shown  in  Figs.  4-6. 

The  data  about  grid  systems  are  summarized  here.  Let  us  put  cf  =  {2kTi/m^)^^'^.  We  divided  the  interval 
0  <  Xi  <  D  into  100  uniform  sections  for  Kn=  1  and  10  and  into  100  nonuniform  sections  (minimum  size 
4  X  10^®D  at  Wi  =  0  and  D\  maximum  size  0.0294Z)  at  Xi  =  D/2)  for  Kn=  0.1.  We  used  uniform  grids  for 
For  Kn=  0.1  and  1,  the  grid  size  is  0.15c^  and  the  range  is  restricted  to  — 6cj  <  <  6cj  {wP  jm^  =  0.5) 

or  <  8.7c^  {mP /mP  =  0.25)  for  A-component  and  to  —SAcf  <  <  8.4c^  {m^ /m^  =  0.5)  or 

—12cf  <  ^1  <  12cf  {m^ /m^  =  0.25)  for  S-component;  for  Kn=  10,  the  grid  size  is  0.106c^  and  the  range 
is  restricted  to  — 5.73c;^  <  <  4.45cj  {m^ /m^  =  0.5)  or  — 7.85cj  <  <  6.58c;^  [m^ /m^  =  0.25)  for 

A-component  and  to  — 7.85cj  <  Ci  <  6.58cj  (jrP  jm^  =  0.5)  or  — 10.82cj  <  <  9.55cj  {vnP /m^  =  0.25)  for 

B-component.  For  r],  we  used  nonuniform  14  grid  points  defined  by  {2kTi /m°‘)^P pijk  {a  =  A,  B)  for  Kn=  0.1 
and  1  and  {2kTi/m'^)^^'^^yk/2  for  Kn=  10,  where  yk  {k  =  1,  ...,  14)  are  the  zeros  of  the  Laguerre  polynomial 
Li4(y)  (see  [10]). 


FIGURE  2.  Profiles  of  the  number  densities  and  n®  and  of  the  temperature  of  the  total  mixture  T  for 
vnP  jrrP  =  0.25,  =  0.5,  and  rp^juav  =  L  See  the  caption  of  Fig.  1. 


FIGURE  3.  Profiles  of  the  number  densities  rP  and  rP  and  of  the  temperature  of  the  total  mixture  T  for 
ra^ /mP  =  0.25,  d///d^  =  0.5,  and  rp^jn^v  =  10.  See  the  caption  of  Fig.  1. 
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(b)  Kn  =  1 


^tl(2kTijm^Yl'^ 

(c)  Kn  =  10 


FIGURE  4.  Velocity  distribution  functions  and  at  two  points  near  the  walls  for  vnP  jm^  =  0.25,  d® /d^  =  0.5, 
and  riavInYv  =  0.1  (cf.  Fig.  1).  (a)  Kn  =  0.1,  (b)  Kn  =  1,  (c)  Kn  =  10.  Here,  —  and  —  indicate  the  result  by  the 
finite-difference  method,  and  •  and  o  that  by  the  DSMC  method.  The  at  X\j D  =  0.095  and  0.905  are  shown  in  (a) 
and  (b),  while  F°‘  at  X\jD  =  0.09  and  0.91  are  shown  in  (c). 
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FIGURE  5.  Velocity  distribution  functions  and  F®  at  two  points  near  the  walls  for  mP  jm^  =  0.25,  d® /d^  = 
and  nav/riav  =  1  (cf.  Fig.  2).  (a)  Kn  =  0.1,  (b)  Kn  =  1,  (c)  Kn  =  10.  See  the  caption  of  Fig.  4. 
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FIGURE  6.  Velocity  distribution  functions  and  at  two  points  near  the  walls  for  mP  jm^  =  0.25,  d® /d^ 
and  =  10  (cf.  Fig.  3).  (a)  Kn  =  0.1,  (b)  Kn  =  1,  (c)  Kn  =  10.  See  the  caption  of  Fig.  4. 


FIGURE  7.  DSMC  computations  with  two  different  numbers  of  simulation  particles  for  vnP  jm^  =  0.25,  d^/d^  =  0.5, 
Kn  =  0.1,  and  =  10.  (a)  Number  density  n'^  and  temperature  T  (cf.  Fig.  3).  (b)  Velocity  distribution 

function  at  Xi/D  =  0.905  (cf.  Fig.  6).  Here,  —  indicates  the  finite-difference  result,  o  the  DSMC  result  with 
(iV^,  iV®)  =  (25,  250),  and  •  that  with  (N^,  iV®)  =  (100,  1000). 


Finally,  we  give  information  about  the  DSMC  computational  system.  We  used  50  (Kn  =  10)  or  100  (Kn  =  0.1, 
1)  uniform  cells  in  the  interval  0  <  Xi  <  D.  Let  AT"  be  the  average  number  of  simulation  particles  per  cell  for 
Q!-component.  Then,  {N^,  iV^)=(1000,  100)  for  =  0.1,  (250,  250)  for  n^y/n^y  =  1,  and  (100,  1000) 

for  n^yln^y  =  10  [(iV^,  iV®)=(2000,  200)  for  n^y/n^y  =  0.1  and  (200,  2000)  for  n^y/n^y  =  10  in  the  case  of 
/m^  =  0.25,  d^/d^  =  0.5,  and  Kn  =  10].  The  average  of  2  x  10^  samples  taken  at  each  50  time  steps 
is  shown  in  Figs.  1-6.  For  small  or  large  n^y/n^y,  the  total  number  of  simulation  particles  increases  because 
sufficient  particles  are  necessary  for  the  component  with  smaller  number  density  (the  same  weight  is  used 
for  both  components  in  the  present  computation).  We  also  carried  out  the  DSMC  computation  with  fewer 
particles,  an  example  of  which  is  shown  in  Fig.  7.  That  is,  the  result  with  {N^,  N^)={25,  250)  of  the  case 
/m^  =  0.25,  d^/d^  =  0.5,  Kn  =  0.1,  and  n^y/n^y  =  10  is  shown  in  the  figure,  together  with  the  result 
with  {N^,  A?’'®)=(100,  1000).  Although  it  is  smooth,  the  profile  of  with  {N^,  N^)={25,  250)  deviates 
recognizably  from  that  by  the  finite-difference  method  [Fig.  7(a)]. 
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